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RECOVERING A POLYHEDRAL OBSTACLE BY A FEW 
BACKSCATTERING MEASUREMENTS 

JINGZHI LI AND HONGYU LIU 


Abstract. We propose an inverse scattering scheme of recovering a polyhedral 
obstacle in n = 2,3, by only a few high-frequency acoustic backscattering 
measurements. The obstacle could be sound-soft or sound-hard. It is shown 
that the modulus of the far-held pattern in the backscattering aperture pos¬ 
sesses a certain local maximum behavior, from which one can determine the 
exterior normal directions of the front sides/faces. Then by using the phaseless 
backscattering data corresponding to a few incident plane waves with suitably 
chosen incident directions, one can determine the exterior unit normal vector 
of each side/face of the obstacle. After the determination of the exterior unit 
normals, the recovery is reduced to a hnite-dimensional problem of determining 
a location point of the obstacle and the distance of each side/face away from the 
location point. For the latter reconstruction, we need make use of the far-held 
data with phases. Numerical experiments are also presented to illustrate the 
ehectiveness of the proposed scheme. 
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1. Introduction 

This work concerns the inverse scattering problem of recovering an impenetra¬ 
ble obstacle by the corresponding acoustic wave detection. The problem has its 
physical origin in radar/sonar, geophical exploration, non-destructive testing and 
medical imaging (cf. (4[[To]). Let D be a bounded Lipschitz domain in n = 2, 3, 
such that M.^\D is connected. D represents an impenetrable obstacle located in 
the space and it is assumed to be nnknown/inaccessible. In order to identify D, 
one sends a time-harmonic detecting plane wave of the form 

«Xx) = e“-«, (1.1) 

which is an entire solntion to the Helmholtz eqnation (—A — = 0 in 

The presence of the obstacle D interrupts the propagation of the plane wave, 
leading to the so-called scattered wave field u^, which exists only in the exterior 
of the obstacle. The total wave field u = satisfies the following Helmholtz 
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system 


(-A- 

lim 

|x|—)-+oo 


iei> = o 

du^ 
d\x\ 


in 


\x\ 


^\D, Bu = 0 on dD, 
= 0. 


( 1 . 2 ) 


In (1.2), Bu := u or Bu := du/dp, with n G S” ^ := {x G 


X 


= 1 } denoting the 


exterior unit normal vector to dD, corresponding to that D is sound-soft or sound- 
hard, respectively. In the following, we set /c = |^| G M+ and d = ^/|^| G 
denoting the wavenumber and incident direction of the plane wave, respectively. 
The PDE system (1.2) is well-understood with u G possessing the 

following asymptotic expansion (cf. E[T^) 


u{x) = -|- 


oik\x\ 


-U 


\x\ 


\x\ 


1 


\x\ 


I '>T- + 1 


I X I —y “h(X), 


(1.3) 


which holds uniformly in x := x/|x| G where x G M” and x 7 ^ 0. m°°(x) is 

known as the far-held pattern and we shall write u°°{x;^, D) = u'^(x;k,d,D) to 
specify its dependence on the observation direction x, wavenumber k and incident 
direction d, as well as the obstacle D. m°°(x) is real-analytic in x, and hence if it 
is known on any open patch of then it is known on the whole sphere by the 

analytic continuation; see |^. 


The inverse problem that we are concerned with in the present paper is to recover 
D by knowledge of u°°, which is known to be nonlinear and ill-posed (cf. |4[[l0]). 
It is noted that the inverse problem is formally posed with a hxed ^ G M” and all 
X G Hence, there is a widespread belief that one can recover D by using 

the far-held pattern corresponding to a single incident plane wave which is 
referred to as a single far-held measurement. However, this still remains to be a 
longstanding problem with very limited progresses in the literature. If the obstacle 
is of small size; roughly speaking, smaller than half of the detecting wavelength, the 
unique recovery result was established in [^. If the obstacle is extremely “rough” 
in the sense that its boundary is nowhere analytic, the unique recovery result was 



obstacles. Another challenging issue in the study of inverse scattering problems is 
the recovery by phaseless data, say the modulus of the far-held pattern, |u°°(x)|. 
The only recent result we are aware of is , where the unique determination of a 
scattering potential by the phaseless far-held measurements was established. 


In this paper, we develop a novel scheme of recovering a polyhedral obstacle by 
using only a few high-frequency far-held measurements. The obstacle could be 
sound-soft or sound-hard. The crux is the observation that the modulus of the 
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far-field pattern in the backscattering apertnre possesses a certain local maximnm 
behavior, from which one can determine the exterior normal directions of the front 
sides/faces. Then by using the modulus of the backscattering data correspond¬ 
ing a few incident plane waves with suitably chosen incident directions, one can 
determine the exterior unit normal vector of each side/face of the obstacle. Af¬ 
ter the determination of the exterior unit normals, the recovery is reduced to a 
finite-dimensional problem of determining a location point of the obstacle and the 
distance of each side/face away from the location point. For the latter recon¬ 
struction, the far-field data with phases would be used. Our study is based on 
the high-frequency asymptotics, namely the Kirchhoff or the physical optics ap¬ 
proximation. However, our numerical experiments show that the high-frequency 
requirement could be relaxed to a certain extent. Moreover, in order to simplify the 
discussion, we focus on convex polyhedral obstacles in the present study. However, 
through our theoretical arguments, it can be expected that the method developed 
would work for non-convex obstacles, but under certain geometrical constraints. 
We focus on developing the novel inverse scattering scheme for convex obstacles in 
the present study and leave the technical and tedious derivation of the geometrical 
conditions for non-convex obstacles for a forthcoming work. 

The rest of the paper is organized as follows. In Section 2, we consider the physical 
optics approximation on the high-frequency scattering from an admissible poly¬ 
hedral obstacle, and derive the local maximum behavior of the modulus of the 
corresponding far-field pattern. In Section 3, we present the recovery scheme in 
detail. Section 4 is devoted to numerical experiments to validate the applicability 
and effectiveness of the proposed method and the paper is concluded in Section 5 
with some discussion. 


2. Physical optics approximation 

Throughout the present section, we let k E M+ and d G be fixed. Let D be a 
convex polygon in or a convex polyhedron in such that 

m 

dD = \JC„ (21) 

i=i 

where each Cj represents an open side/face of dD, and shall be referred to as 
a cell in what follows. In the sequel, D is referred to as a polyhedral obstacle. 
Let u{x) G X G dD denote the exterior unit normal vector to dD, and we 

set 

Uj := u{x) when x E Cj, j = 1,2,... ,m. (2.2) 

Clearly, I'j is a constant unit vector. 
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Define 

dD^ := {x G dD] p(x) • (i > 0} and dD~ := {x G 9D; p(x) • d < 0} 

to be, respectively, the back-face and front-face of dD with respect to the incident 
direction d. 


Let ho, hi and h 2 be fixed a priori positive constants. It is further assumed 
that 


(i) k ■ diam(D) ^ 1; 

{ii) min diam^n-i (Cj) > ho; 

l<j<m 

(Hi) hi < min < h 2 for y G dD. 


(2.3) 

(2.4) 

(2.5) 


Roughly speaking, (2.5) implies that the obstacle should not be very “round” or 
“sharp”, and a generic condition which can guarantee this assumption is that the 
angle between any two adjacent cells is bounded below and above by certain con¬ 
stants (depending on the obstacle). Assumption (i) means that we are considering 
the scattering in the high-frequency regime. If a polyhedral obstacle D satisfies 
the above three assumptions, then it is called an admissible obstacle with respect 
to the incident plane wave 


Denote 


■.= {xeS^-^; x-d>0} and := {x e x ■ d < 0} 

the forward-scattering and backscattering apertures, respectively. 

Let Cj C dD~ be a front-cell of dD, and Vj G denote its unit normal vector 
pointing to the exterior of D. Define Xj G satisfying 

{d - Xj) II Uj 

to be the critical observation direction with respect to d and Pj; see Fig. 1 for a 
2D illustration. We note that one clearly has Xj G It is directly calculated 

that the critical direction is given by 

Xj = d — 2{d • Pj)Pj. ( 2 . 6 ) 


On the other hand, for the subsequent use, we note that by using (2.6) and the 
fact that d- Uj < 0, one has by innerly producting both sides of (2.6) with d 


d-Pj = - 


1 — Xj ■ d 


(2.7) 


Hence, by combining (2.6) and (2.7), one further has 


= 


Xj — d 


^/2(l — Xj ■ d) 


( 2 . 8 ) 
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Figure 1. 2D illustration of the relation between the incident 
direction d, the exterior unit normal vector Pj and the critical ob¬ 
servation direction Xj. 



Figure 2. (a) Plot of the real part of the scattered wave 

in the vicinity of a sound-soft triangular obstacle, and (b) polar 
graph of the square power of the phaseless far-field data in 

polar coordinates corresponding to the triangle due to an incident 
plane wave with d = (—■s/2/2, v^/2). The black arrow indi¬ 
cates both the incident and the forward-scattering directions, while 
the two red arrows represents the two critical observation directions. 
The green arrows indicates the exterior directions normal to the sides 
of the triangular obstacle. 


Next, we present the major result motivating the recovery scheme that we are 
going to develop in the next section. Before that, we give a numerical example 
by plotting the real part of the scattered wave field and the square power of the 
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modulus of the associated far-held pattern corresponding to a sound-soft triangle 
due to an incident plane wave impinging from southeast to northwest; see Fig. 
for the illustration. 

Clearly, one can see a certain local maximum behavior from Fig. [^b). The obser¬ 
vation direction (black arrow) associated with the largest magnitude of phaseless 
data points to the incident/forward-scattering direction, which gives no informa¬ 
tion of the obstacle. While the other two critical observation directions (red arrows) 
associated with the second and third largest magnitude of phaseless data provide 
profound information of unknown obstacle, which signify the major rehection an¬ 
gles due to the physical optics approximation in the high frequency scattering. 
The exterior unit normal direction can thus be determined by connecting the ar¬ 
row heads of incident and critical observation directions by yellow lines. Those 
unit vector (green arrows) parallel to the yellow lines are the desired normal di¬ 
rections of those sides of the triangular obstacle in the backscattering aperture as 
shown in Fig. |^b). We shall present a rigorous mathematical justification for this 
local maximum behavior, which also forms the basis of our subsequent recovery 
scheme. 


We first give some discussion on the Kirchhoff or the physical optics approximation 
(cf. 11, 16,^), which plays a key role in our recovery scheme. Considering 


the present study, it says that for the scattering from a convex polyhedral obstacle 
due to a high-frequency plane wave, the wave field near the boundary of the 
obstacle is composed of two parts: the contribution of the incident and reflected 
waves where they are present, and the contribution of the diffraction due to the 
corners and/or the edges of the obstacle. The Kirchhoff or the physical optics 
approximation takes the first contribution as the total wave field near the boundary 
of the obstacle. However, we would like to remark the study on rigorously justifying 
such approximation is still not fully understood and we refer to for an excellent 
account on the existing progress in the literature. For our study, we assume that 
the physical optics approximation holds true. 


Throughout the rest of the paper, we let a cell C be parameterized as 

(t, u) = I, X e C, (2.9) 

where p G is the unit normal vector to C pointing to the exterior of the 

obstacle, and I > 0 denotes the distance between the origin and the line/plane 
containing C. Let denote the affine cell of C defined by 

{x,p) = 0, x&C^. (2-10) 

In the sequel, we let H^o denote the Euclidean reflection with respect to C^. Now, 
we consider the scattering locally near a boundary cell, say Cj due to a plane wave 

ik(x-xo)-(n^od) ikxn-d 1 ^ ^ 

Vj{x) := e i ■ e ° , 1 < j < m, (2.11) 
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where Xq G Cj is a fixed point, and x G It is easily verified that Vj{x) 

satisfies {A + k‘^)vj = 0. Moreover, one can further verify that 


uUx) — Vj(x) = 0 and ^ q ^ 

ou 


( 2 . 12 ) 


Vj(x) is the reflected wave fleld of «*(x) with respect to the cell Cj. Using the 
physical optics approximation, we have 

Lemma 2.1. Let D be an admissible polyhedral obstacle with respect to the plane 
wave field as described in (2.1)-(2.5). Using the physical optics ap¬ 

proximation, one has 


du 

if D is sound-soft; and 


2du 


du (^)> 

0 , 


X £ Cj C dD , 1 < j < m, 
X G Cji C dD~^, 1 < / < m, 


U[X) 


2n*(x), X £ Cj C dD , 1 < j < 


m, 


0 , 


X e Cj' C dD+, 1 < f < m, 


(2.13) 


(2.14) 


if D is sound-hard. 


Proof. According to our earlier discussion on the physical optics approximation, 
one takes 

^ X eCj C dD~, (2.15) 


dnj 


if D is sound-soft. Then by using ( 2.11[ ), one directly verifles (2.13) for x G dD~, 
whereas if Cj' C dD^, it is not illuminated, where one then takes du/dvjj 0. If 
D is sound-hard, then one takes 

u{x) ~ u{x) + v^(x), X E Cj C dD~, (2.16) 

which readily verifles (|2.14[). □ 


It is emphasized again that Theorem |2.1| is mainly based on physical observation, 
and its rigorous justification still largely remains open in the literature. Never¬ 
theless, our subsequent study on the inverse scattering problem, along with the 
corresponding numerical experiments, validates such approximation as well. 


In the sequel, we let 

gik\x-y\ 

^{x,y) = 


47r|a: — y\ ’ 


n = 3; 


^H^^\k\x-y\), n = 2; x-fiy, 


where denotes the zeroth-order Hankel function of the first kind. $ is the 
fundamental solution to —A — k^. The following lemma shall be needed and its 
proof can be found in |^. 
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Lemma 2 . 2 . For the scattering of a plane wave field r* in (1.1) from an obstacle 
D, we have 


u{x] D, R*) = R*(x) + 


IdD 


f d^(x,g) 
I du{y) 


and 


u°°{x; D, R*) = 7 (n, k) 


.JdD 


f)ii 1 _ 

u(y) - ^x,y)-^(y)> ds{y), i e R”\£l, 

(2.17) 

(2.18) 


^ -ikx-y p! 


where the dimensional parameter 7 is given by 


7(n, k) = — when n = 3; , 

471 


when n = 2. 


(2.19) 


Lemma 2.3. Using the physical optics approximation in Lemma \Kl\ for the scat¬ 
tering of a plane wave in S from an admissible polyhedral obstaele D, one 
has that if D is sound-soft 


«°°(x) ~—27(77,/c) j 

Ja 

whereas if D is sound-hard 


Q^ikyd 

dv{y) 


ds{y), 


( 2 . 20 ) 


u 


r p^-ikx-y 

fix) ^ 2^{n,k) / ^ dsfiy). 

JdD- ^^[y) 


( 2 . 21 ) 


Proof. This is a straightforward consequence of (|2.18|) in Lemma |2.2| and, (pT3|) 
and (2.14) in Lemma 


2.1 


□ 


By using Lemma |2.3[ we are now in a position to present the local maximum 
behavior of \uoo{x; D,k,d)\. 

Theorem 2 . 1 . Let D be an admissible sound-soft or sound-hard polyhedral obsta¬ 
cle with respect to the incident plane wave e*^*'*^ as described earlier. Suppose that 
Cj C dD~ is a front cell of the obstacle, and Uj is the unit normal vector to Cj 
pointing to the exterior of D, 1 < j < m. Let Xj G be the critical observation 
direction with respect to d and Vj. Under the physical optics approximation, Xj is 
a local maximum point of | 7 (°°(x; D, 
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Proof. We first consider the case that D is an admissible sonnd-soft polyhedral 
obstacle. By (2.20) in Lemma 2.3, we have 

r Q^ikyd 


u°°{x; D, k, d) ~ —2'y{n, k) 


^-ikx-y 


'oD- dp{y) 


'y ds{y) 


= — 27 (n, k) / iku{y) ■ ds{y) 


JdD- 

— 2j{n,k) ■ik 


ic, 


Uj ■ d e*^ 2 /-(d-£) ^ f ^ ^ikyid-x) 

JdD-\Cj 


= 7 (n, k) ■ [Ii{x) + hix)], 

where 7 (n, k) := —2^{n, k) ■ ik and 

h{x) := [ Uj-d ds{y), hix) := 

Jci 


( 2 . 22 ) 


'dD-\Cj 


v{y) ■ d ds{y). 


We next analyze the behavior of Ia{x), a = 1, 2, in a small neighborhood of Xj on 
S”“^, say Tj, 1 < j < m. Since 


{d - Xj) II Pj, 


(2.23) 


by assnmption (2.5), we see that there exists a non-asymptotic constant cq E 
such that 


\T{y) • (d - Xj)\ > Co, ye dD \Cj, 


(2.24) 


where T{y) represents a unit tangent vector on dD. Hence, by (2.24) and 
one has by direct calculations that 

, -Cl, X E Tj. 

k^-^ho ^ 


\h{x)\ 


(2.25) 


(2.26) 


On the other hand, for x G Fj, we have 

h{x)= [ Pj ■ d ds{y) 

Jcj 

=pkyo-{d-x) f ^iky(d-x) ^ 

Jq 

where r/o is any fixed point on Cj. Since 

y ■ {d — Xj) =0 for y E Cj, 

one clearly sees that |/i(x)| achieves its local maximum value at Xj, which in 
combination with ( 2.25[ ) completes the proof of the theorem for the sound-soft 
case. 
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The sound-hard case can be shown by following a similar argnment. By nsing 


(2.21), one has 


u°°(x; D, k,d) K. — 2^{n, k) ■ ik / u(y) • x e 

JdD- 

=^{n,k) • [Ji{x) + J 2 {x)], 


iky-{d—x) 


ds{y) 


(2.27) 


where 

Jl(x) = 


'Ci 


Uj'X e 


iky(d-x) _ 


'dD-\Cj 


iy{y)-xd’^y<‘^-^Usiy). (2.28) 


We consider the behavior of Ja{x), « = 1,2, in a small neighborhood of Xj on 
say Sj, 1 < j < m. Following a similar argument to (2.26), together with 


the fact (2.23), one can see that |Ji(t)| achieves its local maximum value at Xj. 
On the other hand, by a similar argument in deriving (2.25), one has that 

\J 2 {x) \ "Cl for T G Ej, 1 < j < m. (2.29) 

Hence, xj is local maximnm point of \u°°(x] D, 

The proof is complete. 

□ 


Remark 2.1. In Theorem 2.1, we only eonsider the local maximum behavior of 


|u°°(f)| in the baekscattering aperture. In fact, by (2.22), one elearly sees that 
|u°°(f)| aehieves its (global) maximum value at x = d. However, this maximum 
behavior in the forward-scattering aperture gives us no information about the ob¬ 
stacle, whieh is demonstrated and verified in Fig. 


Remark 2.2. It is noted from (2.27) that for a sound-hard obstacle if x F p. 


j’ 


then both Ji and J 2 vanish, and hence such an x should be a local minimum point 
of \u°°(x; D,k,d)\. However, it eannot be guaranteed that sueh an x belongs to 
the baekseattering aperture. Nevertheless, it could be used in combination with the 
critical observation direction as a double indieator. 


3. Recovery scheme 


Based on our study in the previous section, we shall present a recovery scheme 
of identifying an admissible polyhedral obstacle D. To that end, we let A: be a 


sufficiently large wavenumber such that assumption |2.3| is fulfilled. Then, let da, 
a = l,2,...,pbea few properly chosen incident directions. The basic requirement 
is that the union of the front-faces with respect to da, a = 1,2,... ,p should cover 
the whole boundary dD. Furthermore, we shall assume that D is an admissible ob¬ 
stacle with respect to each da- As remarked earlier, this assnmption is generically 
satisfied by an overall not very “round” obstacle. 
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First of all, clearly, by Theorem |2.1| and the discussion following its proof, one can 
recover a family of exterior unit normal vectors. However, it must be emphasized 
that multiple groups of normal directions are actually determined corresponding 
to detecting waves with different incident directions. Within each group, those 
normal directions point roughly to the same direction. Up to this stage, we choose 
the normal direction, which is associated with the critical observation angle with 
the largest magnitude of phaseless data within the group, to be an effective exterior 
unit normal direction. Following this principle, we can recover all the exterior unit 
normal vector uj to each of the cell Cj to dD, j = 1,2,... ,m. In turn, we also 
know the number of cells of the obstacle. 


Next, we proceed to identify each cell, and we let xq E D he a. fixed location point. 
Let each Cj be parameterized as follows 

{x — To, Pj) = lj, X e Cj, (3.1) 

where Ij denotes the distance between the origin and the line/plane containing the 
cell Cj — To := {t — To; t G Cj}. We note that Uj and Ij, j = 1,2,..., m, and Tq 
uniquely determine the convex obstacle D. It is highlighted that the location point 
To can be initially guessed by Scheme I in using single-shot measurement data 
at low frequency. Hence, after the determination of To and I'j, j = 1,2,..., m, the 
corresponding inverse problem reduces to a finite dimensional problem of finding 
Ij, j = 1,2,..., m. But we have richer data set u°°(f; k, da) to that purpose. 


Now, let us take the sound-soft obstacle as an illustration to derive the rest of the 
recovery scheme. Suppose t" is a critical observation direction corresponding to 
Cj with respect to a certain da, 1 < a < p. Then, by (2.22|), we have 


u 


jxj; k, d.y) ^ 7(n, k) u, ■ ds{y ). 


(3.2) 


'C, 


Eq. (3.2) is clearly a finite dimensional nonlinear problem with respect to Ij's, 
which are hidden in the implicit functions defining C^’s. 


By supplementing m particularly chosen normal directions and their respective 
critical observation angles, one can arrive at a nonlinear system of m equations 
with m unknown Ij, j = 1,2,..., m. The cells can be determined all at once by 
solving the finite dimensional problem. It is emphasized that the chosen m non¬ 
linear equations (3.2) hold only approximately, thus we convert it into a nonlinear 
least-squares minimization problem to yield a more stable and accurate identifica¬ 
tion. 


Now we are in a position to present our main reconstruction algorithm of polygonal 
obstacles which is sketched in Fig. 


Recovery Scheme. 
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Step 1. Employ Scheme I proposed in 13 to locate the position xq of the polygonal 
scatterer using an incident detecting wave of low frequency. 


Step 2. Determine the m effective exterior normal vectors Pj, j = 1,2, ...,m 
and their associated critical observation angles from multiple groups of candidate 
normal vectors associated with all the critical observation angles. 


Step 3. Given the distances Ij's from Xq along the direction Pj, j = 1, 2,..., m, we 
can locate all the perpendicular points Pj’s. 

Step 4. Extend the line/plane at P/’s perpendicular to the respective Vj and denote 
those crossing points to be vertices Pj’s. The line/plane of the polygonal scatterer 
can be determined by Cj := KjKj/ where j' = j + 1 if j < m and j' = 1 if 
j = m. 


Step 5. Select the first m critical observation angles Z/ 3 , (3 := /3{a,j) = 1,2,... ,m 
from all the x", a = 1,2,... ,p, j = 1,2,... ,m according to the decreasing order 
of magnitude of local backscattering maxima and convert the corresponding m 


nonlinear equations associated with (3.2) with m unknowns l/s into a least-squares 
minimization problem. It is noted that the integral to the right hand side of 


(3.2) is approximated by the trapezoidal quadrature rule with sufficiently fine step 


size. 


mm 


u'- 


m = E 

13=1 

t (b G2; • • • ) ^m) 


{zf^; k, da) - 7(n, k) / Uj ■ da ds{y) 


iCi 


Step 6. The minimum of the multivariable nonlinear least-squares minimization 
problem is obtained by employing a derivative-free trust region method via a local 
quadratic surrogate model-based search algorithm. Interested readers may refer 
to and the references therein. 


4. Numerical experiments and discussions 


In this section, we present some numerical tests to verify the applicability and 
effectiveness of the proposed recovery scheme in two dimensions. In the sequel, 


the forward equation (1.2) is first solved by using the quadratic finite element 


discretization on a truncated circular domain enclosed by a PML layer. The for¬ 
ward solver is iterated on a sequence of successively refined meshes till the relative 
error of two successive finite element solutions between the two adjacent meshes 
is below 0.1%. Then the scattered data are transformed into the far-field data 
by employing the Kirchhoff integral formula on a closed circle (2D) enclosing the 
scatterer. 
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Kj 


Figure 3. Sketch of the recovery scheme. 


First of all, let’s fix the parameter settings. For the positioning pnrpose , we apply 
Scheme I (see [^) with A: = 1 to detect the initial displacement guess xq of the 
unknown polygonal scatterer. A few detecting waves with the incident directions 
chosen among the set {dj = (cos(j7r/4), sin(j7r/4))}, j = 1, 2,..., 8, are sent off for 
locating the sides of the polygonal scatterer. Then the far-field data are measured 
and collected at 360 equidistant observation angles along the unit circle. The 
far-field data generated on the unit circle are then subjected pointwise to certain 
uniform random noise. The uniform random noise in magnitude as well as in 
direction is added according to the following formula, 

u°° = u°° + Sri\u°°\exp{iTrr 2 ), (4-1) 

where ri and r 2 are two uniform random numbers, both ranging from -1 to 1, and 
5 represents the noise level. 

For polygonal obstacles, we shall test sound-soft and sound-hard, noise-free and 
noisy cases, respectively. In the noise-free case, it is well-known that the far-field 
data is analytic and thus very smooth. As a consequence, the local maximum 
behavior of phaseless far-field data is clear from its polar graph as shown in Fig. 
when there exists no noise. While in the noisy case, we always add to the exact far- 
field data a uniform noise of 5% and use it as the noisy measurement data, which 
is inevitable from the practical point of view. But the local maximum behavior of 
phaseless far-field data might be corrupted by the ups and downs of random noise 
which cause fictitious and/or more local mamima than usual. The cure out of the 
dilemma is to add a preprocessing step to filter the raw noisy data. In our tests, a 
fourier filtering stage is applied to the noisy far-field data in advance. The filtered 
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measurement data are thus used as in the recovery scheme, which gives us better 
reconstructions than using raw data by our experience. 



(a) 


(b) 



(c) (d) 

Figure 4. Plot of in polar coordinates corresponding 

to a sound-soft triangle due to an incident plane wave with 

k = 6% and d = d 2 ,d 4 ^,dQ,d 8 from (a) to (d). The selected critical 
observation angles Zi, £2 and £3 are highlighted by red arrows. 


Example 1. A triangle. 


The obstacle is chosen to be a triangle with three vertices displaced at (1,0), 
(2.5, —0.5) and (2.5,1), respectively. In this test, we send off four detecting waves 
from north, east, south and east with d = d 2 , ^ 4 , d^, dg, respectively. 

Firstly, we test a sound-soft triangular obstacle. The initial position xq of the 
polygonal scatterer is detected by Scheme I (see (^) and obtained by taking the 
position with maximum indicator function value. The initial guess of the location 
point is found to be xq = (2.136, 0.217), denoted by a red star in Fig. 
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We plot the square power of the phaseless far-held data in Fig. It can be seen 
from Fig. 1^ that the phaseless data display signihcant maxima along the forward¬ 
scattering (or incident direction) directions within the shadow region of the ob¬ 
stacle in all four plots. Except the forward-scattering directions, we hnd several 
other directions with local maximum behavior in the polar plots and determine 
the three critical observation angles with the largest magnitnde of phaseless data, 
which is indicated by the red arrows in Fig. |^a), (c) and (d), respectively. With 
the respective set of incident and critical observation directions obtained in hand, 
one can determine the three exterior nnit normal directions of the sides of the 
polygonal obstacle. 

The rest of the work is reduced to be a three-dimensional nonlinear least-squares 
minimization problem. The final reconstruction results are shown in Fig. which 
is quite satisfactory in both noise-free and noisy cases. It can be observed that 
we can obtain better reconstrnction by using detecting waves with relatively high 
wave number for this sound-soft polygonal scatterer. This can be explained in 
the following way. The detecting wave with high freqnency gives a more focused 
reflected beam and thus yields a better approximation of the physical optics. The 
higher the frequency of the detecting wave, the more accurate determination of 
the outward normal directions and thus we can obtain better reconstruction plots. 
Moreover, we see that the proposed recovery scheme is tolerable to relatively high 
level of noise and performs robust in the noisy case. 


Next, we keep the experimental settings unchanged except replacing the obstacle 
by a sound-hard triangnlar scatterer. As before, the location point is detected to 
be xq = (1.9307, 0.1412) nsing Scheme I in 13 , denoted by a red star in Fig.[^ As 
shown in Fig.|^we plot the square power of the phaseless far-field data and indicate 
the first three critical observation angles within the backscattering aperture. The 
final reconstruction results are shown in Fig. which performs as good as in the 
sound-soft case. The effectiveness of the recovery scheme can be explained in a 
similar way as in the sound-soft case. 


Example 2. A convex hexagon. 


In the second example, the obstacle is chosen to be a sound-soft hexagon with six 
vertices displaced at (4,2.5), (3,3) (1,2), (0.5,0) (2,-1) and (4.5,—0.5), respec¬ 
tively. The location point is detected to be Xq = (2.582, 0.759) using Scheme I 
This example is much more challenging since there are multiple facets to 


13 


m 


be determined. 


It is pointed out that the hexagonal obstacle has six sides and we only send off 
four detecting waves along incident directions d = di, d^, d^, d-j. In this case, the 
unknown nnmber of sides is larger than that of incident directions. But from Fig. 
one can reveal snfficient critical observation angles from those polar plots. We 
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(a) 


(b) 




(d) 

Figure 5. Reconstruction of the triangular sound-soft scatterer 
with (a) k = Gtv with 5% noise, (b) k = Qtt without noise, (c) 
k = IOtt with 5% noise and (d) k = IOtt without noise. 

identify six critical observation angles with relatively large local maximum values 
as indicated in red arrows in Fig. The final reconstruction results are obtained 
by solving a six-dimensional nonlinear problem and are shown in Fig.|^ It is again 
observed that the reconstruction performs better with higher frequency detecting 
waves. With the combining effect of increasing sides and noise, our recovery scheme 
still yield quite reasonable approximation to the original hexagon. 

5. Concluding remarks 


In this paper, we develop an inverse scattering scheme of recovering a sound-hard or 
sound-soft polyhedral obstacle by only a few backscattering far-field measurements. 
It has been a very challenging issue in the literature on recovering an obstacle by 
minimum measurement data. We believe that the results in this work make some 






















RECOVERING POLYHEDRAL OBSTACLES 


17 



(a) 



(c) 




Figure 6. Plot of \u°°{x)\'^ in polar coordinates corresponding 
to a sound-hard triangle due to an incident plane wave with 

k = 6% and d = d 2 ,d 4 ^,dQ,d 8 from (a) to (d). The selected critical 
observation directions Zi, Z 2 and Z 3 are highlighted by red arrows. 


important contribution to this challenging issue. The proposed scheme proceeds 
with two steps. First, one uses the local maximum behavior of the modulus of the 
backscattering far-held data to determine the exterior normal direction of each of 
the side/face of the obstacle. Then, one can solve a small-scale hnite dimensional 
algebraic problem to completely recover the obstacle. In order to justify the local 
maximum behavior of the modulus of the far-held data, we made essential use 
of the high-frequency asymptotics of the acoustic scattering. Our method can 
be extended to recovering non-convex obstacles, as well as to the electromagnetic 
scattering problems, which we shall report in the forthcoming work. 
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(a) 




(d) 


Figure 7. Reconstruction of the triangular sound-hard scatterer 
with (a) k = Gtv with 5% noise, (b) k = Qtt without noise, (c) k = lOvr 
with 5% noise, and (d) k = IOtt without noise. 
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Figure 8. Plot of in polar coordinates corresponding to a 

convex sound-soft hexagon due to an incident plane wave with 
k = 6 % and d = di,d 3 , 65,67 from (a) to (d). The selected critical 
observation directions Zi, £ 2 , £3, £4, £5 and £5 are highlighted by red 


arrows. 
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(C) (d) 

Figure 9. Reconstruction of the sound-soft hexagonal obstacle 
with (a) k = Gtv with 5% noise, (b) k = Qtt without noise, (c) 
k = IOtt with 5% noise, and (d) k = IOtt without noise. 
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